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We study the quasinormal modes of massless charged fermions in a Reissner- 

^ Nordstrom-anti-de Sitter black hole spacetime. In the probe limit, we find that the 

£>^j imaginary part of quasinormal frequency will become positive when the temperature 

of the black hole is below a critical value. This indicates an instability of the black 

hole occurs and a phase transition happens. In the AdS/CFT correspondence, this 

transition can be viewed as a superconducting phase transition and the bulk fermion 

is regarded as the order parameter. When the coupling of the fermions and the 

background electric field becomes stronger, the critical temperature of the phase 

transition becomes higher. If the interaction between the fermion and the electric 
5_i \ 
(3JT)! field can be ignored, namely in the case of a neutral fermion, the imaginary part 
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of the quasinormal modes is always negative, which indicates that the black hole is 
stable and no phase transition occurs. 



I. INTRODUCTION 

Over the past years the holographic models of superconductors have attracted a lot of 
attentions since the works |l|, |2[ • For reviews see [3j . In a simple model with the Einstein- 
Maxwell-complex scalar field system with a negative cosmological constant, the charged 
scalar field will condensate when the temperature of the anti-de Sitter (AdS) black hole 
is below a critical value. Above the critical temperature, the system has the Reissner- 
Norstrom (RN)-AdS black hole solution with a trivial scalar field, while below the critical 
temperature, a hairy black hole solution is more stable with a nontrivial scalar field. This 
indicates that a phase transition happens between the RN-AdS black hole and a hairy black 
hole when the temperature of the black hole arrives at the critical value. According to the 
AdS/CFT correspondence [4H6J], the phase transition of the AdS black hole can be mapped to 
a superconducting phase transition on the boundary of the AdS space |2| , the condensation 
of charged scalar field around the black hole corresponds to a condensation of the charged 
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operators in the boundary field theory. 

It is well known that a neutral scalar field perturbation of an asymptotically AdS space- 
time is stable if the mass of the scalar field satisfies the Breitenlohner-Freedman (BF) 
bound [7|. However, the perturbation of charged scalar field in an asymptotically AdS 
spacetime can cause the background to become unstable |lLl8MlO||. In the holographic super- 
conductor models, the presence of the instability of the perturbation of the charged scalar 
field just indicates the occurrence of the condensation in the boundary field theory. 

To reveal the stability of a black hole, a useful method is to study the quasinormal modes 



of some perturbations in the black hole background, for a review see JjJ . Recently it was 



found that quasinormal modes of charged scalar field in some AdS black hole backgrounds 



can be used to disclose the relation to the superconductor transition [12|. The occurrence of 
the unstable modes was observed in consistent with the superconducting phase transition. 
Besides the charged scalar field, it is of great interest to examine whether some other 
fields can experience the condensation in the AdS black hole background, for example the 
charged fermionic field, its condensation might be related to the model of color supercon- 
ductor in QCD theory. The perturbation of a charged fermionic field in an asymptotically 



flat spacetime was found always with decay modes [13(. The fermionic field will either be 
absorbed by the horizon or go to the spatial infinity and cannot condensate near the black 
hole. It is well known that the AdS spacetime has different asymptotical behavior from that 
of the asymptotically flat spacetime and this difference plays crucial role on the dynamics of 
the inner geometry. For example, the late-time tail of the perturbations with a power law 
form in an asymptotically flat spacetime gives way to the exponential form in an asymptot- 



ically AdS spacetime [14|, Il5|. In this work we will concentrate on the study of the charged 
fermonic field perturbation in the background of an AdS black hole. 

We employ a simple RN-AdS spacetime geometry while leaving fermionic fields as probes. 
We find that when the temperature is below a certain critical value, the imaginary part of the 
quasinormal modes will become positive, which indicates the occurrence of the instability of 
the black hole. This shows that at the critical temperature there appears a phase transition, 
the RN-AdS black hole background becomes unstable and a new stable hairy black hole 
solution is expected to emerge. We observe that the critical temperature grows with the 
increase of the charge coupling constant. This means that the stronger the fermions couple 
to the Maxwell field, the easier the instability occurs. Our result is consistent with the 
holographic superconducting phase transition discussed in [3j] 12[, however in our result the 



order parameter is a charged fermion instead of the charged scalar field. In this sense our 
model might be regarded as a simple model of color superconductor in QCD theory. 

This paper is organized as follows. In Sec. Hi] we briefly introduce the RN-AdS geometry 
and the action of the charged fermions in this background. Under an explicit representation, 
the Dirac equation is splitted in Sec. [TTT]to a concise form. In Sec. HVl we use the Horowitz- 
Hubeny approach to Fourier expand the Dirac equations for the convenience to study the 
quasinormal modes of the fermionic perturbation numerically. We give the numerical results 
and analysis of the quasinormal modes in Sec. |V] We conclude our paper in Sec. IVIl 



II. THE BACKGROUND GEOMETRY AND ACTION 

For a general d- dimensional RN-AdS spacetime, the metric is 

ds 2 = -f(r)dt 2 + ^ + r 2 dtf d _ 2 , k , (1) 

^ere f(r) = k + g + J^ - (^~ 3 . (2) 

L is the AdS radius, Q is the charge of the black hole and ro is related to the black hole 
mass M via 

d _ 3 _ 16nG d M 
r ° ~(d-2)A d ^ {S) 

where A d ~2 = 27r^~/T(^p) is the area of a unit (d — 2)-sphere, G d is the Newton gravi- 
tational constant in the ci-dimensional spacetime. dfi d _ 2 k in (pQ) is the metric of constant 
curvature. If k = 0, it is the metric of a flat Euclidean space R d_2 ; if k > 0, it is the line 
element of a (d — 2)-sphere with radius -4=; and if k < 0, it is the metric of a hyperbolic 
plane with radius of curvature -7=r- Without loss of generality, one can take A; = 0, and ±1. 
We will focus on k = case in this paper for simplicity and also for the relation to a 
superconductor in a plane. For the cases with fc ^ 0, similar results can be obtained. In 
particular, we will consider a 4-dimensional RN-AdS spacetime and let G d = 1 

dv 

ds 2 = -f(r)dt 2 + -— + r 2 {dx 2 + dy 2 ), (4) 

The temperature of the black hole now is 

H 4vrL 2 16vrr3' W; 

where r + is the horizon radius of the black hole. 
The spin connection is defined as: 

u hhc = e M d c el + e M e{T d fc , (6) 

where e^ is the tetrad *, T a bc is the Christoffel connection and u^ c = —^~ hhc - The nonvanishing 
spin connections uj^ c for the metric (j4j) are: 

bJfit = ^ift = nl ) tdxrx ^fxx v J i ^y?y = ^?yy = y J • \ • ) 

where a prime denotes the derivative with respect to r. 



1 The un-hatted letter b is the index of the background spacetime while the hatted letter a denotes the 
index of the tangent space. 



We consider the Einstein-Maxwell-fermion system with the action 

St = S g + S m , (8) 



where 



1 /"*. ^(- 6 



s 9 = —J d*xV^[n- T2 ), (9) 



S m = M i ' d A xyf=g~ (-^F ab F ab + i(VT a (D a - iqA a )9 - mW)j 



(10) 



Ki is the 4- dimensional gravitational constant, 1Z is Ricci scalar, M is a total coefficient of 
the action of matter, q is the coupling constant between the fermion field and Maxwell field. 
T a is the Dirac gamma matrix and 

T b = 4P\ (11) 

The Dirac equation for the fermion is 

T a (D a -iqA a )m -mm = 0. (12) 

In the action (JSj), clearly we have a simple RN-AdS black hole solution with the electric 
potential A t and a vanishing fermion \l/: 

A t = Q(--—), 9 = 0. (13) 

r r + 

In the following, we will work in the probe limit which means that the fermionic field does 
not backreact on the metric and Maxwell field. 

III. THE SPLIT OF THE DIRAC EQUATION 

We can write the wave function \I/(r, x M ) into the momentum space 

V(r, Xll ) =il>(r)e~ iut+it *. (14) 

where x^ denotes the coordinates in the boundary, x^ = (t,x,y), while k = (k x ,k y ) and 
x = (x,y). Under this transformation, we can write the Dirac equation (|T2"|) into 

Jrd r i> - -^rfy + ^-^ + 7H7? + — ]rV - W a A a + m)v = o. (15) 



where k ■ T x = k x T x + k y Y^. We can choose the Dirac gamma matrices as 16 



r -(°- 'y- °:> ™ 



where / is the 2x2 unit matrix, a 1 is the Pauli matrix, explicitly, 

•*-(;;)• ••-(!?)• '-(;y- ^ 

It is easy to decompose the fermion filed \& into \1/ + and \E'_, which are the eigenvectors of 
T 5 , i.e. 

$ = ( Jj; + V p ± ^ = ±*±, p± = 1 ± r 5 , r 5 = iT t r B r ,, r r . (is) 

Under this representation, the Dirac equation (TT5l) can be decomposed into 

1 f \fl i i 

(y/fd r + -4= + ^)a f ^_ + -(A; • a)^- + -^(w + gA)^- - ™A+ = 0, (19) 

(yffd r + ]4f + —V^ + -$■ ^+ -4f(" + qA t )ip + - m^_ = 0. (20) 
4V7 r r v7 

It is easy to see from eqs.( lT9|) and ( |20l) that if w — >■ —a; and g — > —q, there is a permutation 
symmetry ip + -H- ?/;_ in ( fl9l) and ( |20l) . In the following, we will take m = 0, because the 
chiral modes of ip will decouple in this case, which can be easily seen from eqs. ( 119]) and 
(EQD. 

For the massless fermions, we will focus on the modes of ip + , because the modes of ijj- 
can be obtained accordingly if we change the sign of u and q. Furthermore, for simplicity, 



we can set k y — because of the symmetry of (x, y)-plane [17|. We rewrite ip+ as 

V> + = r- 1 /- 174 ^. (21) 

And then eq. ( 1201) can be further simplified into 



v7, 

p- • r 



a f 4>' --(u + qA t -^±k x a £ )^ = 0. (22) 



Here ib is a 2-component fermion wavefunction, we can decompose it as ib = ( , ] under 
the a r matrix. Under this representation, eq. ( |22l) can be further decomposed into 

ib[ - ~(w + gA)^i + ~^k x ib 2 = 0, (23) 

^2 + t(w + gA)^ 2 - ^7=Mi = 0. (24) 

/ ry/f 



Like in the case of eqs. f fl9|) and (120]) . if tu — > — ui, q — > — q and /c x — )■ —k x , there is also a 
permutation symmetry tb\ o ^2 between eqs. (123]) and (124)) . Therefore, we can only focus 
on the modes ib\ by extracting ib 2 from eq.( |23|) and then inserting it into eq.(J21 



IV. THE QUASINORMAL MODES 

To calculate the quasinormal modes of the charged fermion, we introduce the tortoise 
coordinates r* by 

dn = j. (25) 

Now, the black hole horizon is located at r* — > — oo and the spatial infinity is at a finite 
value of r*. 

The ingoing boundary conditions at the horizon requires the wavefunction to have the 
form 

^(r) = e - iwr *u(r). (26) 

For simplicity of the numerical calculation, we define x = 1/r. Thus, the black hole horizon 
is at x + = l/r + while the infinite boundary r — > oo locates at x — 0. And the equation of 
u is given as 

S( x )u{x) + -^U(x) + . Ui - X \ u(x) = 0, (27) 

x — x + [x — x + ) z 

where the overdot denotes the derivative with respect to x and 



S(x) --- (-L 2 Q 2 x 3 x 3 + + Ax 2 + + Axx+ + Ax 2 ) 2 , (28) 

0). 

(29) 



T(x) = - (L 2 Q 2 x 3 x 3 + - 4x 2 + - 4xx+ - Ax 2 ) (4L 2 Q 2 x 3 x 3 + + 16iL 2 ux 3 + - 3x 2 (L 2 Q 2 x% + 4)) 



U(x) = 2L 2 x 3 + {2uj [-8L 2 qQx A + - UL 2 Q 2 x 3 x% + 8L 2 qQxx 3 + + 3ix 2 [L 2 Q 2 x\ + 4)] 
- {x - x+) [2 (L 2 Q 2 x 3 x 3 + - Ax\ - Axx + - Ax 2 ) k 2 x + iqQ (-8iL 2 qQx\ 
+2L 2 Q 2 x 3 x\ + 8:4 + 8x {iL 2 qQx\ + x+) - x 2 (3L 2 Q 2 xl + A))] } . (30) 

Although the complex number i occurs in the potential U(x) here, it was argued in 18M22I 



and references therein that one still can have correct quasinormal modes in the cases with 
such a complex potential if some proper boundary conditions are imposed. 



Next we employ the Horowitz-Hubeny method [14|] to calculate the quasinormal modes of 
this charged fermion. S(x),T(x) and U(x) can be Fourier expanded near x + to a finite term, 
i.e., S(x) = Yl=o s n(x - x + ) n ,T(x) = Y? n =o f n( x ~ X +T and u (x) = Yl=Q u n(x- x + ) n . 
We also expand the solution u(x) around x + as 

oo 

u{x) = (x-x + ) a ^a n (x-x + ) n . (31) 

n=0 

At the leading order, the equation has two solutions a = 1/2 and a = , 2 -t 2 o^ a u - ^ ^ s 



12-L 2 Q 2 x 4 , 



8iL 2 x 



obvious that the solution a = 12 _l 2 q^ 4 u gi yes an outgoing mode at the horizon. As a 



result, we choose a = 1/2 in the expansion (131]) . Substituting eq. (l3~Tj) into the eq. (!27|) . we 



can obtain a recursion relation 22 



n— 1 

B" ^[(fe + a)(A; + a - l)s n _ fc + (fc + a)t n _ fc + u„_ fe ]a fc , (32) 

*" fc=0 



where P n = {n + a){n + a — 1)sq + {n + a)t + u$. 

In the numerical calculations, we set the initial data ao = 1 whose scaling will not affect 
the final quasinormal modes u because the equation of ip is linear. In practice, we can solve 
the expansion up to a large number n = N, then compare the results got from N and the 
results from n > N, if the error of the two results is of the order 10~ 5 , we will believe that 
the modes un is acceptable. 

V. NUMERICAL RESULTS 

The boundary condition on the spatial infinity is that the fermion field vanishes there, 
which means 



N 
n=0 



(0) = ^a n (0-x + r +Q = 0. (33) 



This equation gives the eigenvalues of u. The quasinormal modes can be decomposed into 
real and imaginary parts: 

cu = Re(w) + ilm(u). (34) 

If Im(o;) < 0, the perturbations always decay exponentially and then the background black 
hole is stable; if Im(o;) > 0, however, the modes grow exponentially, which means the 
perturbations will make the background black hole unstable. This instability suggests a 
phase transition of black holes [8|, [12] . 

Following pj, in the calculation we set r + — Q — 1. In this case, changing L corresponds 
to the change of the temperature of the black hole. Furthermore, we set the x-direction 
momentum k x — 1. For the calculation precision we Fourier expand the equation up to 
N = 300. 

To find the marginally stable perturbation where \I/ depends only on r and is infinites- 
imally small, taking u = in fl27|) . we will have an equation on q and L from eq. fl27|) . 
Requiring both the real and imaginary parts of the resulting equation on q and L to vanish, 
we have two constraints on the relation between q and L, whose numerical results are shown 
in the left panel of FigJU It indicates that only q and L at the crossing points of the two con- 
straint curves can give the marginally stable solution with u = 0. The three crossing points 
(q,L) shown in the left panel of Fig fl] are (1.546, 1.951), (2.821, 1.915), and (4.187,1.887), 
respectively. When q increases, we see that L at the crossing points decreases. According 
to the temperature formula ()5]) of the black hole, the temperature increases as L decreases. 
Thus the stronger the fermion field couples to the Maxwell field, the easier the instability 




Real Part 
Imaginary Part 



FIG. 1: (Left Panel) Eq. (|27|) with u) = gives the two constraints on q and L: the real part (blue) 
and imaginary part (red dashed). (Right Panel) The critical temperature Th c versus q. 



occurs. This is physically reasonable. We show this relation of the critical temperature and 
the coupling q in the right panel of FigfTJ 

The marginally stable perturbation cannot guarantee the existence of the instability but 
just highly suggestive [l|. In examining the stability of the perturbation, one can release the 
condition u = 0, then in principle there does not exist any constraint relation between q and 
L. In this paper we are interested in the quasinormal modes of the fermion perturbation and 
expect to see whether a positive imaginary part of quasinormal frequency which indicates 
the instability of the fermion perturbation will appear. Figj2] shows the frequencies of the 
quasinormal modes in the case with a fixed q = 1.546 and different L. It can be seen that 
when L > L c = 1.951 the imaginary part of the quasinormal frequency becomes positive, 
which means that the black hole is unstable under the perturbation when the temperature 
is low enough. L = L c corresponds to the critical temperature of the black hole Th c , when 
T > Thc the RN-AdS black hole is stable under the fermonic perturbation, while the RN-AdS 
black hole becomes unstable when T < Th c and the fermonic hair is expected to condensate 
and attach to the black hole. The RN-AdS black hole will give way to a fermion haired 
black hole. 

In Table [J we list the quasinormal modes shown in Fig.2. One can see that the imaginary 
part of the first node becomes positive when L = 1.96 and L = 2, while it is negative if 
L = 1.94. This manifests that the perturbation becomes unstable when L changes from 1.94 
to L = 1.96. The critical value is L = 1.951, where u = 0. 

From the left panel of Fig. [1] we can see that when q = 0, there are no critical values for L 
which makes u = 0. This implies that for the neutral fermion perturbation, the black hole 
is always stable and the perturbation always decays exponentially. This is indeed observed 
in the quasinormal modes calculation as shown in Table [Til where we see that the imaginary 
part of the quasinormal modes of the fermion perturbation is always negative in the case of 
q = 0. 

Fig. [3] shows the behavior of the quasinormal modes with respect to the black hole tern- 
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FIG. 2: The real and imaginary parts of the modes when changing L while q is fixed as q = 1.546. 



TABLE I: Real and imaginary parts of the quasinormal modes shown in Fig[2j 



Node 


1st 


2nd 


3rd 


4th 


5th 


6th 


7th 


8th 


9th 


10th 


L = 2 


Re(cj) 


0.0402 


0.4437 


1.2797 


2.1952 


3.1451 


4.1073 


5.0763 


6.0094 


6.9258 


7.8056 


Im(w) 


0.2311 


-0.0581 


-0.5066 


-0.9778 


-1.5873 


-2.3222 


-3.1731 


-4.1324 


-5.1931 


-6.3485 


L = 1.96 


Re(u>) 


0.0407 


0.6464 


1.7115 


2.8367 


3.9965 


5.1664 


6.3280 


7.4656 


8.5658 




Im(w) 


0.1136 


-0.2701 


-0.7109 


-1.3015 


-2.0632 


-2.9793 


-4.0383 


-5.2309 


-6.5482 


L = 1.94 


Re(w) 


0.0665 


1.5676 


2.8389 


4.1304 


5.4067 


6.6680 


7.9068 








Im(w) 


-0.1366 


-0.6918 


-1.0962 


-1.8305 


-2.7393 


-3.7939 


-4.9829 



perature. We can see from Fig. [3] that in the high temperature regime, both the real part 
and imaginary part of the quasinormal modes have a linear relation to the temperature of 
the black hole, which is similar to the result given in 22[, where the quasinormal modes of a 
neutral massless fermion are numerically calculated in a RN-AdS black hole with a spherical 
horizon (namely the case with k — 1 in ([I])). The linear behavior observed for the fermonic 
perturbation agrees to the finding for scalar perturbation [l4J, |l5j. By numerical fitting, we 
find that the quasinormal frequencies behave like 



Re(u) « 8.4193T H , Im(u) » -6.3981T H , 



(35) 



TABLE II: Real and imaginary parts of the quasinormal modes in the case of q 
temperature shown in Fig. [3l 



with different 



T H 


0.0537 


0.0733 


0.1019 


0.2188 


0.3531 


0.6432 


1.4721 


2.6326 


5.9484 


10.5904 


Re(w) 


0.8253 


0.9856 


1.2209 


2.2134 


3.3206 


5.7491 


12.6862 


22.3972 


50.1418 


88.9839 


Im(w) 


-0.5725 


-0.6992 


-0.8804 


-1.5724 


-2.4708 


-4.3155 


-9.5919 


-16.9810 


-38.0943 


-67.6533 



10 



where T# is the black hole temperature given in (jSJ). In the high temperature regime, 
Ref. 22j gives a relation (see eq. (8) in 22] ) for the case of Schwarzschild-AdS black hole 
with a spherical horizon 



Re{u) = 8.367T, Im{u) = -6.371T. 



(36) 



We can see from ( 1351) and ( |36l) that their behaviors are quite similar to each other, except 
the differences in the numerical factors caused by the value of the nonzero black hole charge 
Q we considered here. 

In the low temperature regime, on the other hand, one can see that both the real and 
imaginary part obviously deviate from the linear behavior. This result is also consistent 
with the one of scalar fields in [14|, Il5| . 




FIG. 3: Quasinormal modes of the first node in the case of q = 0. The real part versus high 
temperature (Up-Left); the imaginary part versus high temperature (Up- Right); the real part versus 
low temperature (Down-Left); the imaginary part versus low temperature (Down- Right). 



VI. CONCLUDING REMARKS 



In this paper we considered the Einstein-Maxwell-Fermion system with a negative cos- 
mo logical constant. Such a system has been intensively studied recently in the holographic 



fermion liquid model [16|, llll I23K25J . We studied the stability of the RN-AdS black hole 



11 

under the perturbation of the charged fermion by numerically calculating the quasinormal 
modes of the perturbation in the probe limit. It is found that when the temperature of the 
black hole is below a critical value, the imaginary part of the quasinormal modes will become 
positive, which means that the perturbation of the fermion will grow exponentially which 
makes the black hole unstable; while when the temperature of the black hole is above the 
critical value, the perturbation will decay exponentially and the black hole is stable. Our re- 
sult indicates that a phase transition occurs in the system from the high temperature phase, 
which is described by the RN-AdS black hole, to a low temperature phase, which should be 
described by some stable new black hole solution. Like the case with charged scalar field, we 
expect that the stable new black hole should be a one with nontrivial fermion condensation. 
When the coupling constant q increases, the critical temperature of the black hole becomes 
higher. This means that the stronger the coupling is, the much easier for the phase transition 
to occur. In the case of q = 0, namely for a neutral fermion, we found that the imaginary 
part of the quasinormal modes is always negative. This implies that the black hole is 
stable under the perturbation of the uncharged fermion and it will not cause any phase 
transition of the black hole. The behavior of the quasinormal modes agrees with that 
for the scalar perturbation [l4j, |l5] and in consistent with fermonic perturbation 22[ for 



spherical background. In high temperature regime, both the real and imaginary parts of the 
quasinormal modes are linearly proportional to the temperature of the black hole. In the 
AdS/CFT correspondence, l/Im(o;) is the time scale for the system approaching to thermal 
equilibrium of the boundary field theory. This means if we perturb the thermal field on the 
boundary, the time for it to approach a thermal equilibrium is proportional to the inverse 



of the temperature [14J. However the linearity breaks in the low temperature regime. 

In the probe limit, we found the marginal stable mode exists with u = of the charged 
fermion. This strongly indicates there should exist a stable charged black hole with fermion 
hair in AdS space. It would be of great interest to consider the backreaction of the fermion 
and to find such a black hole solution in our system, which is currently under investigation. 
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